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NON-EQUILIBRIUM  HOT-CARRIER  DIFFUSION  PHENOMENON  IN  SEMICONDUCTORS 
I.  A  THEORETICAL  NON-MARKOVI AN  APPROACH + 


J.  Zismermann* ,  P.  Lugli  and  D.K.  Ferry 

Colorado  State  University,  Fort  Collins,  CO  80523 ,  £/.S.d. 

* Permanent  address  :  C.B.S,  and  Greco  Microondes,  University  de  Lille  1,  France 

Res xmik  -  Le  p rob 1 feme  de  la  diffusion  de  port curs  chauds  dans  les  a eaal con¬ 
duct  eur  a  est  fetudife  A  l'alde  d ' une  Equation  de  Langevln  Retardfee  (RLE) 
appllqude  au  cas  de  la  rdponse  transltolre  dynanlque  (TDR)  A  un  champ  dlectrique 
statlonnalre  homogfene.  Une  fonctlon  de  correlation  des  fluctuations  de  vitesse 
non- 8 tat ionnalre  est  ddflnle  et  est  rellde  A  un  coefficient  de  diffusion 
dependant  du  temps.  Ceci  est  appllqud  au  Slllclum  type-N  et  les  auantltds 
interessantes  sent  dtudides  en  fonctlon  de  l'espace  et/ou  du  temps.  Le 
problfeme  de  diffusion  non-statlonnalre  est  partlcullferement  Important  dans 
les  components  A  canaux  ultra-courts  oh  rdponae  transltolre  dynamlque  et  sur- 
vltesse  se  nan lies tent. 

Abstract  -  The  problem  of  hot  carrier  diffusion  In  semiconductors  Is  studied 
with  a  non-Markovlan  Retarded  Langevln  Equation  (RLE)  applied  to  the  case  of 
carriers  in  the  transient  dynamic  response  (TDR)  to  a  steady  homogeneous  elec¬ 
tric  field.  A  non-st at ionary  velocity  fluctuation/correlation  function  is 
defined  and  related  to  a  time  dependent  diffusion  coefficient.  This  Is  applied 
to  n-type  silicon  and  the  parameters  of  Interest  are  studied  as  a  function  of 
space  and/or  time.  The  problem  of  non- stationary  diffusion  Is  particularly 
Important  in  very  short  channel  devices  in  which  TDR  and  velocity  overshoot 
occur. 


1.  Introduction.-  In  recent  years ,  much  Interest  has  centered  upon  the  transient 
dynamic  response  of  electrons  as  It  Impacts  carrier  transport  through  small  spatial 
regions  of  high  electric  field.  With  recent  improvements  In  technological  fabrica¬ 
tion  of  very-short-chaimel  devices,  this  problem  has  become  not  only  of  theoretical 
Interest  but  of  practical  Interest  as  well.  For  Instance,  In  the  pinch-off  region 
of  a  short-gate  field-effect  transistor,  the  carriers  Injected  at  the  source  move 
by  a  combination  of  drift  and  diffusion  In  a  very  high  electric  field.  Then  the 
transit  time  of  the  carriers  under  the  gate  can  be  shorter  than,  or  of  the  same 
order  of  magnitude  as,  the  time  needed  to  establish  s  steady-state  high-fleld  dis¬ 
tribution  function.  In  fact,  a  condition  for  this  to  occur  Is  t**at  the  transit¬ 
time  In  the  high- field  region  be  comparable  to  the  momentum  relaxation  time  thus 
causing  the  velocity  to  Increase,  but  much  shorter  than  the  energy  relaxation  time. 
Thus,  on  average  the  carriers  may  transit  through  a  considerable  portion  of  the 
hlgh-fleld  region  with  almost  their  low  field  mobility  even  though  the  applied 
field  corresponds  to  the  saturated  velocity  range  (1].  This  Is  true  not  only  for 
the  first-order  moment  of  the  distribution  function  of  the  carriers  (drift  velocity)* 
but  also  Is  true  for  higher  order  moments,  and  especially  for  diffusion  (related  to 
the  second  moment).  The  diffusion  coefficient  Is  one  of  the  most  important 
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parameters  required  In  modeling  semiconductor  devices,  it  Is  not  only  necessary  for 
evaluating  operating  characteristics  and  frequency  characteristics  but  It  provides 
also  fundamental  characterization  of  velocity  fluctuations  In  the  system  and  their 
contribution  to  noise  In  the  device  12 J.  Diffusion  actually  is  a  process  depending 
upon  velocity  correlation  and  the  relationship  between  diffusion  and  drift,  as 
expressed  by  the  Einstein  relation,  is  a  steady-state  relation  [3].  Indeed,  tract¬ 
able  results  for  the  steady-state  hot  electron  problem  have  only  recently  been 
achieved  [4-6].  The  problem  in  the  transient  region  is  complicated  by  the  fact  that 
the  random-walk  equations  governing  transient  diffusion  do  not  reduce  to  normal 
Pick's  lav  behavior  on  time  scales  comparable  to  the  relaxation  process,  a  result  of 
the  general  non -Markovian  and  non-stat ionary  nature  of  transport  on  these  time 
scales.  In  this  paper,  we  address  some  of  these  problems  with  the  help  of  a  Retard¬ 
ed  Langevin  Equation  (RLE)  in  order  to  approach  the  random  walk  of  the  carriers. 
Furthermore,  we  define  a  non-stat ionary  two-time  correlation  function  for  the  ve¬ 
locity  fluctuations  which  can  be  related  to  a  transient  diffusion  coefficient.  The 
formal  solution  of  the  RLE  allows  us  to  derive  a  general  expression  for  the  correla¬ 
tion  function.  Then,  we  deal  with  the  diffusion  coefficient  itself  and  show  in 
particular  that  the  transient  diffusion  coefficient  Is  related  to  the  time  deriva¬ 
tive  of  the  mean-square  displacement  of  the  carriers.  However,  it  is  found  that  in 
the  limit  of  long  times,  statlonarlty  and  the  normal  equations  for  correlation  and 
diffusion  are  recovered. 

2.  The  Retarded  Langevin  Equation  Approach.-  We  consider  an  ensemble  of  carriers 
initially  at  equilibrium  with  the  crystal  lattice.  The  ensemble  is  represented  by 
a  Maxwel  1-Bolt smann  distribution  in  v-epace  (<v>  -  0  and  <v2>  -  3k^T^/m  with  the 
lattice  temperature).  At  a  certain  time,  raf erred  as  t  ■  o,  we  apply  a  microscopi¬ 
cally  homogeneous  and  steady  electric  field  sfcose  amplitude  corresponds  to  hot  car¬ 
riers.  These  conditions  give  rise  to  a  TDR  regime  In  which  the  system  relaxes  toward 
a  non-equilibrium,  steady-state  and  often  exhibits  a  velocity  overshoot. 

The  motion  of  the  particles  la  governed  by  a  Retarded  Langevin  Equation  for  the 

velocity  of  the  carriers ,  and  la  writtan  aa  [7,8] 
t 

■  ft  “  v(u)du  +  R(t)  +  qIoh(t)  ,  (1) 

o 

where  R(t)  la  a  random  force  symbolising  the  random  (non-regular)  part  of  the  col¬ 
lisions  of  the  carriers  with  the  lattice  (no  carrier-carrier  interaction  la  consid¬ 
ered  bars),  Eo  is  the  external  field,  and  h(t)  la  the  Hesvyside  function.  y(t)  la 
the  memory  function  of  the  system  and  Is  related  to  the  correlation  function  of  the 
total  fore*  applied  to  the  system.  For  Instance  in  the  case  of  a  stationary  regime, 
wa  would  have  y(t)  *  <l(s)t(t)>/s2<^(o)>,  In  absence  of  external  forces  (31. 

Equation  (1)  Is  a  nan  Markovian  form  of  the  Langevin  aquation,  since  the  rata 
of  change  of  the  velocity  at  t  not  only  dap  ids  on  the  velocity  at  that  time  but 
alas  depends  on  all  paat  time.  Furthar,  It  la  a  noa  stationary  aquation*  since  the 


lower  bound  in  the  Integral  refers  to  that  tine  where  the  disturbing  field  was 
applied.  It  Is  generally  admitted  by  now  that  only  an  equation  such  as  (1)  can  de¬ 
scribe  very-fast  processes  [3,7,8] ,  and  in  particular  the  TDR  refine  can  be  de¬ 
scribed  in  this  way.  In  (1),  we  have  assumed  a  parabolic  energy  band  [9). 

We  nay  easily  solve  (1)  using  Laplace  transforms.  We  Introduce  a  function 
X(t)  defined  by  Its  Laplace  transform,  as 

X(.)  -  (s  +  V(»))_1  .  (2) 

Then,  coming  back  to  time  domain,  we  find 

t  t 

l  XMin  » j  t(t-a)  X(u)du  ,  (3) 

o  o 

which  Is  a  general  expression  of  the  evolution  of  the  velocity  of  each  carrier  under 
the  Influence  of  the  external  field  and  of  the  collisions.  He  get  X(t)  by  aver¬ 
aging  (3)  over  the  ensemble  (we  assume  <H(t)>  *  0).  The  result  Is 

• 

where  v^(t)  Is  the  ensemble  drift  velocity  of  the  carriers.  Therefore  X(t)  repre¬ 
sents  in  fact  the  macroscopic  acceleration  of  the  ensemble.  However  X(t)  can  be 
given  another  definition.  We  define  a  non-st at lonary  correlation  function  for  the 
velocity  fluctuations  as 


♦Av(t,’°  *  *Av(t,t,)  "  <v<t>v<t,)>  V°Td(t,) 


Multiplying  both  sides  of  (1)  by  v(o)»  averaging  over  the  ensemble  (note  that 
<R(t)v(o)>  -  0),  Laplace  transforming  and  making  use  of  (2)  we  obtain  after  re 
transforming 


♦Av(o,t)  *  <°>>  *<*>  •  <6> 

X(t)  Is  the  reduced  non-st  at  lonary  correlation  function  calculated  et  t*  •  o.  Com¬ 
paring  this  equality  with  (4)  gives 
t 

VO  -  — af2 —  [  ♦.<•,**)*'  .  (7) 

m<v2(o)>  l 

This  relationship  existing  between  the  first  and  a  second  moment  of  the  velocities 
of  the  carriers  le  a  direct  Intrinsic  property  of  the  MLR  need  to  describe  the  Dt 
regime.  Whether  this  equality  le  met  is  erne  of  the  goals  of  the  meat  paper  [10], 
hat  (7)  is  a  statement  of  the  familiar  Ubo  formula  [11]  fowmd  la  eqwlllbrlwm 
statistical  mechanic?. 

To  develop  a  complete  expression  for  the  correlation  function  defined  In  (5)» 
we  need  to  know  the  correlation  function  of  the  random  fores  t(t)  which  appears  whoa 
we  put  (3)  into  (5).  If  we  assume  that  the  collisions  occur  loot  set  ensrrsolj  Is  time 
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<R(t)B(t,)>  -  2*IR(t,t,)«(lt-t,l) 


(8) 


Relaxation  of  this  condition  does  not  significantly  affect  the  results  found  In  the 
present  work.  A  consequence  of  (8)  Is  that  only  I^(t)  *  ^(t.t'-t)  la  of  Interest 


In  the  current  context  and  this  latter  function  can  be  obtained  fro*  the  time  evolu¬ 
tion  of  the 

put  in  the  for*  (with  0  >.  0) 


energy  of  the  ensemble.  Then  the  correlation  function  4.  can  be 

BW 


♦iv^o^o  +  e>  ■  <v2(o)>X(t„)X(to  +  e)  +  |  IR1(to-y)X(y)X(y  4  9)dy  . 


and  In  particular  we  obtain  for  the  mean-energy : 


<e(t)>  -  \  »  vd2(t)  +  \  m  <v2(o)>  X2 (t)  +  -  |  ^(t-y)  X2(y)dy  . 


(9) 


(10) 


He  can  sun  up  this  part  by  saying  that  If,  In  the  current  conditions  of  the  TDt 
regime,  the  evolutions  of  the  drift  velocity  and  the  mesa* energy  are  known,  then  the 
fluctuations  of  the  system  are  completely  specified  through  (9)  and  (10)  together 
with  (4). 

3.  The  Diffusion  Coefficient.-  Using  expression  (5)  for  the  non-stat lonary  correla¬ 
tion  function  of  velocity  fluctuations,  we  can  define  a  non-stat lonary  or  time-de¬ 
pendent  diffusion  coefficient,  which  la  a  generalisation  of  the  definition  given  In 

a  stationary  regime  [3,6]: 
t 


D(t) 


•  | 


t)dt* 


(U) 


Another  way  to  define  the  diffusion  la  through  the  spreading  of  a  packet  of 
carriers  drifting  under  the  Influence  of  the  external  force  and  spreading  due  to  the 
fluctuations  of  the  velocities  of  the  carriers.  This  spreading  Is  characterised  by 
the  mean-square  displacement  of  the  carriers  starting  from  a  known  Initial  distribu¬ 
tion  (a  Dirac- function  In  space,  for  Instance).  It  la  easy  to  see  that 

t  t 
t 

(12) 


<Ax2(t)>  -  <(x(t)-<x(t)>)2>  -  |  } 


o  o 

From  the  definition  (12),  It  la  straightforward  to  show  that 
t 


<4X2(t)>  -  ^♦Av(t',t)dt’  -  »<t) 


(13) 


This  shows  that  the  diffusion  coefficient  defined  in  (11)  Is  related  to  the  time- 
derivative  of  the  sverage-squsre  dlsplacsmsnt .  The  usual  definition  of  the  diffu¬ 
sion  la  a  stsady-state  Is  la  fact  a  limiting  case  of  relation  (13).  Hhan  the  time 
t*  at  which  the  integration  of  4^(t9,t)  begins  la  grsater  than  the  tlms  t#  needed 
for  the  syntax  to  roach  etatloaarlty,  ^  becomes  sn  even  function  of  t-tf  only, 
then 

1  dt  m  \  ♦lT(n)*«  -  !>g<t)  (»> 
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As  t  •,  I)  (t)  tends  to  a  constant  finite  Unit  I)  and 
a  a 

11*  <AK2(t>>  -  2t  (15) 

t  • 

Therefore  (11)  defines  a  general  diffusion  coefficient  which  Is  valid  in  both  sta¬ 
tionary  and  non-statlonary  regimes,  and  as  we  did  above  for  we  can  derive  an 
2 

expression  for  <Ax  (t)>: 

2  t 

<A*2(t)>  -  -f— j  <r2(e)>  »d2(t)  +  f  I^t-u)  »<J2(u)do  .  (16) 

q  o  o 

The  expression  of  D(t)  follows  lasttd lately  using  (13)* 

A.  Application  and_ Discussion In  the  case  of  n-type  S1<111>,  the  carriers  behave 
as  If  the  energy  band  was  unique  with  an  Isotropic  effective  ness  *  ■  There,  the 
concepts  described  above  are  easily  applied. 

1)  The  oscillatory  nature  of  the  velocity  overshoot  strongly  suggests  that 
the  memory  function  y(t)  (and  f(t)v  as  well)  is  an  exponential.  So  we  specify  y(t) 
via  the  ansats  [12] 

Y(t)  -  [Yo  +  »(1  -  .-rt)Jh(t)  ,  (17) 

1/r  la  the  time  needed  for  the  system  to  reach  statloosrlty  and  as  such  Is  equiva¬ 
lent  to  the  energy  relaxation  time* 

11)  To  have  an  Insight  Into  the  time  dependence  of  ^(t),  we  can  Mke  use  of 
(10)  and  the  results  given  by  a  Monte  Carlo  method  [10]  for  <t(t)>  .  It  Is  thus 
possible  to  show  that 


hi(t>  -  **>  +  (v-  Va  -  •  •*> 


r  i  r 

0 


Is  a  good  approximation  for  the  correlation  function  of  the  random  force  E(t). 
Is  given  by  the  static  diffusion  coefficient 


«(•  ♦  ror 

“d  «  Ij^.  The  value  of  1^  does  not  affect  elgnlf leant ly  the  results  since  at 
short  times  the  second  term  In  the  KHS  of  (16)  is  an  order  of  m^ltude  imllr  than 
the  flret  term. 

Taking  (17)  and  (18)  Into  account,  and  !><t)  are  easily  derived  analy- 

tlcally  and  computed.  Me  present  here  eoma  competed  results  obtained  with  ■  -  50 
kV/cn.  the  non-statlonary  correlatioa  function  la  displayed  la  Tig.  1,  and  its 
evolution  Is  studlsd  for  different  Initial  tinea  t0  from  Q  to  0.5  pa  (the  latter 
corresponds  to  the  steady-state) .  The  fact  that  the  correlation  function  la  steeper 
whan  t0  la  longer  can  be  explained  by  the  Increase  of  the  scattering  frequency  as  a 
function  of  time  corresponding  to  the  heating  of  the  carriers.  In  other  words,  at 
short  tinea  the  correlation  function  la  dominated  by  the  initial  distribution  of  the 
velocities  of  the  carriers,  while  at  longer  times,  when  the  steedy-etete  la 
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approached •  the  system  has  been  cca^  lately  randomised  by  the  collisions  and  the 
correlation  function  Is  dominated  by  the  random  force  (i.e.  I^(t)). 

The  mean-square  displacement  and  diffusion  coefficient  are  displayed  In  Tig.  2 
as  a  function  of  the  distance  travelled  by  the  carriers*  This  shows  the  spatial 
extension  over  which  the  system  Is  In  a  non-stat  ionary  regime »  *  400  A  In  the 
present  case*  These  results  are  compared  to  what  would  occur  in  case  of  a  pure 
ballistic  regime.  In  fact,  the  calculated  D(x)  departs  from  the  ballistic  trend 
even  at  very  short  times  (and  distances)  meaning  that  no  ballistic  regime  exists  for 

transient  diffusion.  On  the  contrary  at  longer  times  D(x)  goes  through  a  maximum 

and  then  decreases  to  Its  stationary  value.  This  Is,  of  course,  related  to  the 
oscillatory  nature  of  the  correlation  functions  and  these  oscillations  are  essen¬ 
tially  the  consequence  of  the  combination  of  momentum  and  energy  relaxation  In  the 
resolvent  X(s)  of  the  RLE.  In  the  present  example,  this  resolvent  Is  a  rational 
fraction  which  has  two  complex  conjugate  poles. 

Possible  extensions  of  the  work  could  be:  1)  the  application  of  the  present 
technique  to  multivalley  semiconductors;  11)  to  take  into  account  spatial  varia¬ 
tions  of  the  electric  field  and  then  get  a  more  precise  picture  of  what  occurs  in  a 

short-channel  device;  111)*  to  try  a  more  physical  and  less  phenomenological  approach 
using  quantum  statistical  mechanics  f  131  * 

In  aiwmary,  we  have  obtained  here  a  consistent  definition  of  the  diffusion 
coefficient  In  terms  of  the  velocity  auto-correlation  function.  This  definition  Is 
valid  In  the  transient  sad  non- stationary  regime  and  reduces  to  the  normal  expres¬ 
sion  as  steady- st ate  la  approached.  The  application  of  this  to  the  retarded,  non- 
stationary  Langevln  equation  yields  expressions  for  the  velocity  correlation  func¬ 
tion  and  diffusion  coefficient  which  have  excellent  qualitative  agreement  and  satis¬ 
factory  quantitative  agreement  to  results  obtained  by  a  Monte  Carlo  method  flO]. 
U<«nncM 

1.  H.  Browner,  Sol.  State  Electron..  21  (1978)  81. 

2.  J.-P.  Bougler,  "Boise  Diffusion  in  Semiconductors , r'  In  Physics  of  Bon-Llnear 
Transport  in  Semiconductors.  Edited  by  D*  E.  Perry,  J.  R.  Barker,  end  C.  Jaco- 
bonl  (Plasma  Press,  BswTork)  1979,  pp.  415-464. 

3.  R*  Kabo,  "Response,  Relaxation  end  Fluctuation,"  In  Lectures  Botes  In  Physics  3k 
Transport  Phenomena.  Edited  by  J.  Ehlers,  X.  Rapp,  and  H.  A. . Bsidenmuller 
(Springer  Terlag,  Berlin)  1974,  pp.  74-126. 

4.  C.  Jseobool,  C.  Cornell,  C.  Ottevlanl,  and  A.  Alberlgl-Queranta,  Sol.  State 
Electron..  20  (1977)  77. 

5.  E.  Psaquembergue,  J.  Elamarmanm,  A.  Xaaiynski,  ami  E.  Constant,  J.  AppI.  Phra. 

51  (1900)  1065* 

6.  9.X  ferry,  and  J.  R.  Barker,  J,  Aasl.  Pfrva.  *  52  (1901)  818. 

7#  E.  Swemalg,  Ti  ttfc  OBOU  90  (1960)  2527. 

8.  E.  Borl,  Prom.  Thaar.  ftw,,  33  (1905)  423. 

9.  Bn  cam  take  into  see  pant  mom-parabolic  lty  of  tha  anargy  band 


by  Oaal  tag  with 


DUMoaV') 


C7-101 


the  momentum  p(t)  Instead  of  the  velocity  v(t),  and  the  foraula  v(t)  ■  p(t) 

2  -1/2 

<1  +  2<*p  (t)/m  )  /«  for  exaaple.  In  this  case,  the  derivations  would  be 

o  o 

such  more  Involved. 

P.  Lug 11,  J.  Zlnaermann,  and  D.  K.  Perry,  these  proceedings. 

R.  Rubo,  J.  Phva.  Soc.  Japan.  12  (1957)  570.  See  also  H.  J.  Kreuzer,  Won- 
ulllbrlun  Themodvnaaic  a  and  its  Statistical  Poundatlons  (Oxford  On  tv.  Press 


Oxford)  1981. 

YQ  is  the  inverse  mentua  relaxation  tine  at  thernal  equilibrium  and  a  +  Yq  is 
the  inverse  momentum  relaxation  tine  when  non-eqifl'ilriud  iCMr-dtite-if - r-r 

N  if  7  : 

D.  K.  Perry,  Invited  Paper,  these  proceedings.  I  feL/r  •  !  «5 


\ 

Et  SOkVcur4#^  ii> 

300  K 

*\ 

\\ 

\\ 

Transient  Reduced  Correlation 
functions  Calculated  At 


to-0( — )  ;to-0.02ps(- 


-);t  -0.5ps(-.-) 

o  J 


figure  2. 

Mean-Square  Displacement  (right  scale) 
and  Diffusion  Coefficient  (left  scale) 
as  a  Pune t ion  of  Distance.  Dashed 
Curses  Stand  for  a  Pure  Ballistic 
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